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Abstract 

We have computed the local electric field in a body-centered tetragonal 
(BCT) lattice of point dipoles via the Ewald-Kornfeld formulation, in an 
attempt to examine the effects of a structure transformation on the local 
field strength. For the ground state of an electrorheological solid of hard 
spheres, we identified a novel structure transformation from the BCT to the 
face-centered cubic (FCC) lattices by changing the uniaxial lattice constant c 
under the hard sphere constraint. In contrast to the previous results, the local 
field exhibits a non-monotonic transition from BCT to FCC. As c increases 
from the BCT ground state, the local field initially decreases rapidly towards 
the isotropic value at the body-centered cubic lattice, decreases further, reach- 
ing a minimum value and increases, passing through the isotropic value again 
at an intermediate lattice, reaches a maximum value and finally decreases to 
the FCC value. An experimental realization of the structure transformation 
is suggested. Moreover, the change in the local field can lead to a generalized 
Clausius-Mossotti equation for the BCT lattices. 
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I. INTRODUCTION 



When a strong field is applied to a composite medium, the induced change of the medium 
can lead to spectacular behavior, both in electrical transport and in optical response If 
a strong electric field is applied to a suspension of particles where the particles have a large 
electric polarizability, the induced dipole moments of the particles can order the suspended 
particles into a body-centered tetragonal (BCT) lattice 0. This process is known as the 
electrorheological (ER) effect. 

Recently, Tao and coworkers proposed that a structure transformation from the BCT 
ground state to some other lattices can occur when one simultaneously apphes a magnetic 
field perpendicular to the electric field and the polarized particles possess magnetic dipole 
moments. Sheng and coworkers verified the proposal experimentally and observed a 
structure transformation from the BCT to face-centered cubic (FCC) lattices. Motivated 
by these studies, we propose an alternative structure transformation from the BCT to the 
FCC structure, through the application of electric fields only. 

The plan of the paper is as follows. We will adopt the point-dipole approximation p| 
and calculate the dipole lattice sum via the Ewald-Kornfeld formulation in section II. In 
section III, we discuss the effects of a structure change on the local field when the lattice 
constants vary. If we change the uniaxial lattice constant under the hard-sphere condition, 
a series of transformations occur among the BCT ground state, body-centered cubic (BCC), 
intermediate and FCC lattices (see Fig.l). The results will be compared with those of 
a tetragonal lattice. We also make a contact with macroscopic concept and derive the 
Clausius-Mossotti equation. In section IV, we will compute the dipole interaction energy 
as a function of various lattices and discuss a possible structure transformation by the 
application of rotating electric fields. Discussion and conclusion on our results will be given. 
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II. FORMALISM 



In this section, we apply the Ewald-Kornfeld formulation [^,|| to compute the local 
electric field for a BCT lattice of point dipoles. The BCT lattice can be regarded as a 
tetragonal lattice, plus a basis of two point dipoles, one of which is located at a corner and 
the other of which at the body center of the tetragonal unit cell. The tetragonal lattice has 
a lattice constant c = along the z-axis and lattice constants a = b = iq~^^'^ along the 
X and y axes. The volume of the tetragonal unit cell remains Vc = as q varies. In this 
way, the degree of anisotropy of the tetragonal lattice is measured by how q is deviated from 
unity and the uniaxial anisotropic axis is along the z-axis. 

So far, the lattice parameter ^ remains arbitrary. For hard spheres in an ER solid, 
however, the lattice parameter ^ can be determined from the relation: 2a^ + = 16i?^, 
where R is the radius of the spheres. The hard sphere condition requires a > 2R and 
c > 2R. For the BCT ground state, a = b = ^/QR and c = 2R, while for the FCC lattice, 
a = b = 2R and c = y/8R 0. The BCC lattice is characterized by a = 6 = c = ^16/3R. 
Each sphere has a point dipole embedded at its center. 

The lattice vector of the tetragonal lattice is given by 

R = ^{q-^/^l^ + q-^/^my + qnz), (1) 

where l,m,n are integers. Suppose there are point dipoles pj located at rj in a unit cell. 
The local electric field Ej at a particular point dipole at can be expressed as a sum of the 
electric field of all dipoles at r^j-. 

E. = EET.R.-P., (2) 
j R 

where "prime" denotes a restricted summation which excludes j = i when R = and 

J^, = -V,V,--^ (3) 

is the dipole interaction tensor. Eq.(0) can be recast in the Ewald-Kornfeld form [§,0: 
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3 R 

- ^ II ^ ^^P [(P^ ■ ^) exp(iG ■ Yi) ^(pj- ■ G) exp(-2G ■ r^)] 

+ (4) 

where rjRj = jr^ — trjI, r/ is an adjustable parameter and G is a reciprocal lattice vector: 

G = — (g^/^Mx + g^/^t;y + g^^wz). (5) 

The B and C coefficients are given by: 

erfc(r/r) 2r/ 2 2n 

= 1 — + exp -r] r , (6 

C(r) = 5f£5M + (is; + ^) exp(-,V^), (7) 

where erfc(r) is the complementary error function. Thus the dipole lattice sum of Eq.(^ 
becomes a summation over the real lattice vector R as well as the reciprocal lattice vector 
G. Here we have considered an infinite lattice. For finite lattices, one must be careful about 
the effects of different boundary conditions 0. 

We should remark that although a tetragonal lattice is considered, Eq.(^ is applicable to 
arbitrary Bravais lattices. The adjustable parameter r] is chosen so that both the summations 
in the real and reciprocal lattices converge most rapidly. In what follows, we will limit 
ourselves to the BCT cell with two dipoles per tetragonal cell, and the Ewald-Kornfeld 
summation [Eq.(^)] can be carried out. We will consider two cases depending on whether 
the dipole moment is parallel or perpendicular to the uniaxial anisotropic axis. In both 
cases, we will compute the local field as a function of the degree of anisotropy q. 



III. EFFECTS OF STRUCTURE TRANSFORMATION ON THE LOCAL FIELD 



Consider the longitudinal field case: p = pz, i.e., the dipole moments being along the 
uniaxial anisotropic axis. The local field E at the the lattice point R = reduces to 
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and Ex = Ey = 0. In the equation, and Rj are respectively given by: 



J-1 p 
Zj — n , it j 



j - 1 ^ 
R (ax + ay + cz) 



and ^(G) = 1 + exp[i{u + v + w)/TT] is the structure factor. The local field will be computed 
by summing over all integer indices, {j,l,m,n) ^ (1,0,0,0) for the summation in the real 
lattice and {u,v,w) ^ (0,0,0) for that in the reciprocal lattice. Because of the exponential 
factors, we may impose an upper limit to the indices, i.e., all indices ranging from — L to L, 
where L is a positive integer. For g 7^ 1, the regions of summation will be rectangular rather 
than cubic in both the real and reciprocal lattices. The computation has been repeated for 
various degree of anisotropy with q ranging from 0.5 to 2.0. A plateau value for E^ is found 
for each q within a certain range of 77 values: 1 < 77 < 10°'^. For instance, the calculations 
with ri = 10°'^ yield numerical results already accurate up to 16 significant figures, indicating 
that convergence of the local field has indeed been achieved with the upper limit L = 4. 
For larger anisotropy (either g<^lorg^l),a spherical region of summation can help the 
convergence [^. 

For the transverse field case in which the dipole moments are perpendicular to the uni- 
axial anisotropic axis, Eq. (8) can still be applied to evaluate the local field by modifying Gz 
to Gx while taking the gradient along the direction of the dipole, say the x-axis, and obtain 
the expression of the local field. 

The results of the local field strength (normalized to AirP/S) against log^g*? the 
longitudinal and transverse field cases are plotted in Fig. 2(a) and Fig. 2(b) respectively. For 
comparison, the corresponding results for a tetragonal lattice [§] (i.e., in the absence of 
the body centers) are also plotted on the same figure. As q decreases, the local field for 
the longitudinal field case increases rapidly while that for the transverse field case decreases 
rapidly. In both cases, when q deviates from unity, the effect of anisotropy has a pronounced 
effect on the local field strengths. 



Unlike the tetragonal case, the local field of the BCT lattice does not vary much near 
g = 1. When we magnify the scales in Fig. 2(c) and Fig.2(d), we observe a non-monotonic 
behavior as q increases: as c increases from BCT, the local field initially decreases rapidly 
towards the isotropic value at BCC, decreases further, reaching a minimum value and in- 
creases, passing through the isotropic value again at an intermediate lattice, reaches a max- 
imum value and finally decreases to the FCC value. The isotropic value of the intermediate 
lattice is attributed to the symmetry of the dipole interaction tensor. 

Our present theory is of microscopic origin, in the sense that we have computed the 
lattice summation by the Ewald-Kornfeld formulation. We have not invoked any macroscopic 
concepts like the Lorentz cavity field ||9HTI| in the calculations. However, to corroborate with 
these established concepts can lead to a modification of the Clausius-Mossotti equation. 

More precisely, we use the result of the local field to evaluate the effective polarizability 
OefT of the dipole lattice. The total field acting on a dipole is the sum of the applied field Eq 
and the local field due to all other dipoles, hence 

p = a{Eo + (3P), 

where a is the polarizability of an isolated dipole and (3 = E / P is the local field factor. We 
will use (3z and I3xy to denote the local field factors parallel and perpendicular to the uniaxial 
anisotropic axis. Note that I3z = (3xy = ^iii/'i when q = 1. Let P = p/Vc, the above equation 
becomes a self-consistent equation. Solving yields 

The effective dielectric constant ees is given by 1 + inaes/Vc. Thus, 

= ^ (10) 

where /?' = 3/5/47r and = Vc/2. For the BCC lattice, = 1, satisfies the well- 
known Clausius-Mossotti equation. Thus Eq.(lO) represents a generalization of the Clausius- 
Mossotti equation to the BCT lattices. 
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IV. STRUCTURE TRANSFORMATION VIA ROTATING ELECTRIC FIELDS 



The results of the local field strength as well as the generalized Clausius-Mossotti equa- 
tion allow us to compute the dipole interaction energy per particle in the ER solid, similar 
to calculations of Tao and coworkers [0]. If an uniaxial field E = zi?^ is applied, the dipole 
interaction energy per particle is given hj u = — p - E/2e2, where €2 is the dielectric constant 
of the host medium. Since = PzC^p/Vc) and p = aE^, we obtain 

= ~77TT^ ill units of (11) 

As q increases from the BCT ground state to the FCC structures, the volume of the unit 
cell increases initially, reaching a maximum at the BCC structure, then decreases to the 
FCC structure. As the local electric field strength remains almost constant in the range 
1 < g < 1.5, the increase in the magnitude of the dipole energy per particle is attributed to 
the decrease in the volume of the unit cell. Concomitantly, as shown in Fig. 3, the energy per 
particle initially increases from the BCT ground state, reaching a maximum near the BCC 
structure and then decreases all the way towards the FCC structure. The transformation 
involves climbing up an energy barrier beyond which the FCC structure becomes stable, 
which is in contrast with the smooth and monotonic transition proposed in previous work 

m- 

Next we apply a rotating electric field Exy = rEz in the plane perpendicular to the 
uniaxial electric field, where r = E^y/ E^ is the ratio of the rotating to axial field strength. 
The instantaneous electric fields are E^ = ExyCos{ujt + 0) and Ey = ExySin{ujt + 0), where 
u) is the angular velocity of the rotating field and (p is an arbitrary phase angle. In this case, 
the dipole interaction energy per particle is modified to: 

^ = ~iT^^\iTrt^ ^^^^^ °^ — 3- (12) 

(1 + r^)(Vc/a'^) €20^ 

As shown in Fig. 3, when we increase the ratio r, the BCT ground state energy increases 
while the FCC energy remains unchanged, but there is still an energy barrier between the 
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BCT and FCC states for small r. When r > 2.6, the energy barrier disappears and the FCC 
structure is the most stable. Physically one has to break the chains to make the structure 
transformation possible when the strength of the rotating electric field is sufficiently large. 

From the energy consideration, we suggest that the structure transformation be realized 
in experiments by applying a rotating electric field in the plane perpendicular to the uni- 
axial electric field. In this field configuration, both the time average value of the induced 
dipole moment and that of the rotating electric field vanish in the plane. However, the 
instantaneous dipole moment will induce an overall attractive force between the particles in 
the plane perpendicular to the uniaxial field. This is very different from the previous trans- 
formation proposed by Tao et al and experimentally verified by Sheng and coworkers 
1^ by using crossed electric and magnetic fields on microparticles which possess permanent 
magnetic dipole moments. In our case, the field configuration is all electrical; no magnetic 
field and/or magnetic materials need to be used. Alternatively, we may apply the same 
rotating electric field configuration to a magnetorheological fiuid to achieve the structure 
transformation. 

V. DISCUSSION AND CONCLUSION 

Here a few comments on our results are in order. As mentioned in Ref 0, the energy 
difference between the FCC and hexagonal close-packed (HCP) structures is very small, 
and is comparable to the thermal energy, there is a competition between these structures. 
Therefore, we will likely find an FCC-HCP mixed structure. However, it appears that our 
field configuration helps a four-fold symmetry in the plane of the rotating field and the FCC 
structures may be more favorable. Nevertheless, we are awaiting experimental evidence on 
the proposed structure transformation. 

Our calculations have been performed up to the point-dipole approximation, which is 
in the same spirit as done by Tao et al 0. A dipole-induced-dipole (DID) model, which 
takes into account the mutual polarization effect between touching particles, can drastically 
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improve the accuracy towards the fully multipolar calculations [|T2|. Since the DID contribu- 



tion becomes important for small reduced separation a = d/2R < 1.1 [Q, we can simplify 
the calculations by limiting ourselves to touching spheres only. We may omit the nontouch- 
ing spheres (unless the spheres get very close so that a = d/2R < 1.1). Each sphere is in 
contact with 8 neighboring spheres in the BCT lattices and 12 neighboring spheres in the 
FCC lattice and the number of DID images dipoles is therefore finite. The generalization 
thus includes the original point dipoles as well as all the DID images dipoles at well-defined 
positions in the unit cell and the general formula [Eq.(4)] can indeed be used. The incorpo- 
ration of the more accurate DID model into the Ewald-Kornfeld formulation is underway. 
However, we believe that the point-dipole results remain qualitatively correct. 

In conclusion, we have applied the Ewald-Kornfeld formulation to a BCT lattice of point 
dipoles to examine the effects of structure transformation on the local field distribution. 
We have found that the local field exhibits a non-monotonic transition from BCT to FCC. 
Moreover, we showed that the change in the local field can lead to a generalized Clausius- 
Mossotti equation for the BCT lattices. 



ACKNOWLEDGMENTS 

This work was supported by the Research Grants Council of the Hong Kong SAR Govern- 
ment under grant CUHK 4284/OOP. K.W.Y. acknowledges the hospitality received during 
his participation in the Workshop on Soft Matter, hosted by the Institute of Physics of 
the Chinese Academy of Sciences, where the present work was initiated. We acknowledge 
discussion with Jones T. K. Wan on various aspects of the project. 



9 



REFERENCES 

[1] For recent work, see Proceedings of the 5th International Conference on Electrical Trans- 
port and Optical Properties of Inhomogeneous Media, Physica B 279, (2000). 

[2] R. Tao and J. M. Sun, Phys. Rev. Lett. 67, 398 (1991). 

[3] R. Tao and Q. Jiang, Phys. Rev. E 57, 5761 (1998). 

[4] P. Sheng et al, Physica B 279, 168 (2000). 

[5] P. P. Ewald, Ann. Phys. (Leipzig) 64, 253 (1921). 

[6] H. Kornfeld, Z. Phys. 22, 27 (1924); thesis, Goettingen (Germany), 1923. 

[7] M. Allen and D. Tildesley, Computer Simulation of Liquids, (Oxford Science, London, 
1990). 

[8] C. K. Lo, Jones T. K. Wan and K. W. Yu, in Abstracts of the Conference on Computa- 
tional Physics (CCP2000), December 2000; J. Phys.: Condens. Mat. 13, 1315 (2001). 

[9] C. J. F. Bottcher, Theory of Electric Polarization, Vol.1 (Elsevier, Amsterdam, 1973). 

[10] C. J. F. Bottcher and P. Bordewijk, Theory of Electric Polarization, Vol.11 (Elsevier, 
Amsterdam, 1978). 

[11] H. A. Lorentz, The Theory of Electrons, (B. G. Teubner, Leipzig, 1909). 

[12] K. W. Yu and Jones T. K. Wan, in Proceedings of the 8th International Conference on 
Discrete Simulation of Fluid Dynamics, Comput. Phys. Comm. 129, 177 (2000). 



10 



FIGURES 

FIG. 1. A sequence of unit cells of BCT lattices during structure transformation: (a) Tao's 
BCT lattice, (b) BCC lattice, (c) intermediate lattice and (d) FCC lattice. 

FIG. 2. (a) The local field factor plotted against log;^o ^ dipole moments along the uni- 
axial anisotropic axis, (b) Similar to (a), but for dipole moments perpendicular to the uniaxial 
anisotropic axis, (c) and (d): Magnified versions of (a) and (b) respectively to show the different 
lattices during structure transformation. 

FIG. 3. The dipole interaction energy per particle plotted against q for different ratio 
r = Exy/Ez of the rotating to axial field strength. For hard spheres, the accessible regions are 
within the vertical lines depicted by Tao's BCT and FCC. Along Tao's BCT line and from bottom 
to top, r = 0.0, 1.0, 1.6, 2.0, 2.2, 2.4, 2.6, 3.0, 4.0 and oo. It is evident that the FCC structure is 
the most stable for large r. 
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(a) Tao lattice unit cell 



(b) bcc unit cell 




(c) the Intermediate lattice unit cell (d) fee unit cell 

Fig.l/Lo and Yu 
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Fig.2/Lo and Yu 
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Fig.3/Lo and Yu 
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